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NATIONAL AERONAUTICS AND SPACE ADMINISTRATION 


COMPUTER PROGRAM FOR DETERMINING MASS PROPERTIES OF A 
RIGID STRUCTURE 

BY Reid A. Hull, Dr. John L. Gilbert and Phillip J. Klich 


SUMMARY 


A computer program (ND0702) has been developed for the rapid 
computation of the mass properties of complex structural systems. The 
program uses rigid body analyses and permits differences in structural 
material throughout the total system. It is based on the premise that 
complex systems can be adequately described by a combination of basic 
elemental shapes. The following thirteen widely used structural shapes 
were selected for inclusion in the program: 

1. Discrete mass 

2. Cylinder 

3. Truncated cone 

4. Torus 

5. Beam (arbitrary cross section) 

6. Circular rod (arbitrary cross section) 

■/ . bpherical segment 

8. Sphere 

9. Hemisphere 

10. Parallelepiped 

11. Swept Trapezoidal Panel 

12. Symmetric Trapezoidal Panels 

13. Curved Rectangular Panel 

Simple geometric data describing size and location of each element 
and the respective material deniity or weight of each element are the 
only required input data. From this minimum input, the program yields 
system weight, center of gravity, moments of inertia and products of 
inertia with respect to mutually perpendicular axes through the system 
center of gravity. The program also yields mass properties of the 
individual shapes relative to component axes. 

Permanent configuration records and the use of iterative calculations 
to investigate design systems or to determine optimums contribute to the 
cost-effectiveness of the programs use. 


INTRODUCTION 


Determining the mass properties of any rigid structure is a problem 
that at times becomes complex, but one which can easily be dealt with 
utilizing computer solutions. 



For rigid structures the solution of the mass properties requires 
transformation to an axis parallel to the system axis and becomes 
laborious almost to the point of being impractical. 

Any complex structure must be broken down into elements in order to 
exact a solution. The approach selected for the program presented in 
this paper was to automate the input to the point where an element’s 
shape, geometry, density or weight, and three grid points are the only 
requirements. This approach was influenced by the simplicity of com- 
puting the direction cosines (Euler angle relationship) from the given 
three grid points. The program as outlined in this paper performs 
essentially the same process as calculations "by hand" and is extremely 
useful for rigid structures skewed in space. This program also provides 
improved accuracy, time savings, and complete permanent records for a 
mass properties analysis. (This TMX is a verified expansion of the LaRC 
working paper "Computer Program for Determining Mass Properties of a Compo- 
site Body", by Phillip J. Klich and John L. Gilbert dated Oct. 22, 1968.) 


SYMBOLS 


r 


T 

V 

x’, y’, z’ 
X, y, z 

LJ 


moments and products of inertia 
displacement vector for differential mass 
kinetic energy 
linear velocity 

rectangular coordinate component displacement vectors 
rectangular coordinate system displacement vectors 
angular velocity 
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Subscripts 

CO system coordinates to component center of mass 

xxcc“^ 

yyco ) refer to component axis to vhich moments of inertia are calculated 
zzco J 


XXCO*^ 

YYCO > refer to system axis to which the moments of inertia are rotated 
ZZCO J parallel to the system coordinates 


Superscripts 

prime denotes component coordinate system 


INERTIA EQUATIONS 


The mass properties of shapes such as a cylinder, sphere, etc., are 
easily calculated and therefore were selected as the basic component shapes 
^ 0 '*^ handling a system such as a spacecraft structure. Since the component 
shape mass properties are measured with respect to their respective center of 
masSj the e properties have to be transferred to the system center of mass* 
The transformation can be made in two steps: first, the component properties 
are transferred to a system parallel to the system axis, and then transferred 
by the usual parallel axis theorem. The rotational transformation is derived 
by using the principle of kinetic energy. An introductory derivation of the 
moment of inertia, and product of inertia expressions are derived first and 
then transformation from the component to system coordinates is presented. 


Derivation of Inertia Equations 

Using the expressions of kinetic energy of a rigid body, the equations 
of moments of inertia and products of inertia__are. derived. Consider a 
component spinning with an angular velocity o) as shown next 
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The angular kinetic energy can he written as 



where the velocity is expressed as 

V = cjo X r 


Substituting in the kinetic energy expression gives 

T = — / {(jo X r^ • (a? x r^dm 

2 ^ 


where 


and 


to = io;,. 




+ teB„ 


r = ix' + jy* + kz' 


Taking the cross product 


i j k 
'~‘x ' ‘‘V ' 


= i(oyz' - ci-^y') + JCcu^x' - cu^z') + k(ui3jy‘ 


- o^x') 
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and performing the dot product resxilts In 

(tu'z')^ - 2(cu's ')(a)'y') + (cu'y*)^ + (o)’x')^ - 2 (oj'x’ )( o)’z ’ ) 
'y y Z Z Z Z jC 

+ (cii'z')^ + (‘U’y’)^ - 2(a>'y’)(u>'x') + 

' X X X y y 

Which upon substituting into kinetic energy equation gives 

T = i rf^'^Cy'^ + z'^) + ‘u'^(x’^ + z'^) + u)’^(x'^ + y’^) 

2 J \x y ^ 


_ 2cu'u)'x'z' - 2<u'u)’y'z' - 2u)'o)'x'y' 
z X y z*' X y 


dm 


This represents the rotational kinetic energy of one component of a system. 
Recognizing the definitions of moments and products of inertias and selecting 
the component coordinate system as the principal axes , we can write 

comp 2 ' X X y y z z 


This particular selection of coordinates does not affect the final answers 
because kinetic energy is constant with regards to the coordinate orientation; 
however, it does simplify the input data and also reduces computer time. 


Rotating From Component Coordinates to System Coordinates 

At this point we have the kinetic energy about the component axis system. 
In the computer program it is at this level that mass properties are computed 
for the preselected shapes such as the cylinder, sphere, etc. 

It is necessary to resolve the component mass properties to an axis 
systam parallel to the system coordinates. Once parallel to this system axis 
then we can translate to the system center of gravity by the usual parallel 
axis theorem. The derivation is similar to that presented in reference 5. 

In deriving the rotational transformation from the component to the system 
coordinates the expression for kinetic energy is again used. Given a body 
rotating with an angular velocity ve know that its kinetic energy is 
invariant with regard to the coordinate orientation. 
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sys 




dzn 


“-filing the 

Igg - / da, ly^ . /• 

the energy equation becomes 


yz dm^ 


/ 


xy dm 


+ ly^ + + 2Iyj,oyn^ + 2I^^yO^ 
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It is necessary to write the angular velocity of the cooitonent system 
coordinates in terms of the system coordinates. For an arbitrary vector it can 
be written 

a' = a'* = u' . 

comp system 


Performing the dot product gives 


tL!'=l''iu) + i'*Ju) + i'-lsu 

X X y z 


a'^' = j- . io)^ + j' . jo)y + j’ . 


(n'=k'-itu + k'-jcu +k'-toc 


Recognizing the direction cosines results in 
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Writing the component kinetic energy in matrix algebra 

^comp - |(“')^ . 5 [[DC] {“>; [I'l [tDC]{“}J 

where [DC] is the direction cosine matrix shown above. Equating the system 
and component energy results in 

T ^ = T 

system conQ) 

i tl] {“} = i [DC]'^ [I’] [DC] {a>} 


and the resulting system inertial matrix is found to be 
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[I] * [DC]^ [I’] [DC] 


Expanding the matrix equation we get 


^co 

■ScYCO 

■Sczco 

^YXCO 

^YYCO 

^YZCO 

^ZXCO 

^ZTCO 

^ZZCO 


x’x y*x z’x 

I I I 

x*y y'y z^y 

X 'Z y ^z Z ^Z 

The direction cosines are determined 


0 

0 


Z 

x*x 

1 

x'y 

1 

x»z 

I* 

y 

0 


1 

y’x 

1 

y’y 

1 

y »Z 
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I ’ 
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1 

z *x 

1 , 
z y 

1 , 
z ’z^ 


by the method presented in appendix A. 


I* 

X 

0 

0 


Tre-nsf erring :y :^arallel Axis Theorem 

Nov that the moments of inertia are in a system parallel to the system 
coordinates, we now translate by the parallel axis theorem 

^ ^CO 

lyy = - X)^ + - Z)2] 

"zz " ^ZZCO 

^XT " ^XYCO ■ ^0 

^XZ = %ZCO ^ ■ X)(Z^Q - Z)^ 

^YZ = ^YZCO ■*■ “C^co ■ - Z)^ 
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Where I^q^co* ^C0‘ ^ZZCO *=*** compo^nt moments of inertia 

rotated parallel to the system coordinates, and X, Y, Z are the system 
center of mass coordinates and X^^, Y^^, Z^^ are the coordinates of the 
COTtponent center of mass. 

COMPUTER PROGRAM 


This computer program is written in Fortran IV computer language. All 
namoa and descriptions are assigned in the first part of the program. Thirteen 
sections have been written using 13 common shapes usually found in spacecraft. 
The program is directed by the input data which singles out the secWon 
or shape factor desired to be used through the "go to" statement. The 
operation of the program is illustrated in figure 1 with a computer flow 
diagram. 


The input data for each tiem is listed on two data cards. The basic 
input for each item will vary depending on the shape factor used. Each 
shape factor with the necessary data is discussed in the input data 
instructions . 


After the data cards are supplied to the program the following 
operations are performed. The component mass properties are first printed 
with the moments of inertia about the component axis rotated parallel to the 
system coordinates. These mass properties are transferred to the system 
center of gravity and the following are computed: System weight, inertias 

az iz the system center of gra ^it: , inertias about the origin, center of 
gravity of the system and products of inertias of the system. Based on this 
generated information, inertias about the system principal axes and their 
location is subsequently computed. 

A listing of the computer program is found in appendix B. 


Selection of Coordinate Points 

The selection of points "i" and "j" determines the length of the member 
as well as the first three direction cosines. Point "k" is required to 
calculate the other six direction cosines. Shown in figure 2 are the two 
coordinate systems used in this program. Point "i" locates the system 
coordinates (X., Y , Z. ) for the origin of the component axes and point j 
determines the^direction of the "x" axis of the component coordinates. 

In order to determine the directions of the "y" and "z" axes, point "k 
is required. This point can be anywhere in the x-y plane. If it is omitted 
then the program automatically positions the "y" axis parallel to the X-Y 
plane. For a body of revolution point "k" is not required. The following 
figures (2(a) and 2(b)) describe points 1, j, and k. 


The main deck Is referred to as the computer program without the necessary 
data. It is always necessary to have a 789 card following the main deck and 
a 789 card, then a 6789 card following the data. The 6789 card separates one 
program from another. A description of these cards and their formats follows 
figure 2. 9 
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LOCATION OF POINTS "1" AND *j" 
Figure 2(a) 



LOCATION OF POINTS "k" 
Figure 2(b) 
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1st Data Card 


Item IVb. 

15 format 

Columns 1 through 5 

Description 

2A9 format 
Columns U through 21 

Shape 

12 format 
Columns 22 and 25 

Weight or 
density 

format 

Coliunns 2k through 52 

A 

F8.5 format 
Columns 33 through ko 

B 

PS. 3 format 
Columns kl through kQ 

C 

ER.3 format 
Columns U 9 through 56 

D 

F 8.3 format 
Columns 57 through 6h 

F 

F3.3 format 
Columns 65 through 72 


It is to be noted that the input data variables A, B, C, D, and F can be 
geometric dimensions, cross-sectional areas, area moments of inertia, and 
mass moment of 


2nd Data Card 


XI 

P 8.3 format 
Columns 1 through 8 

YI 

F 8.3 format 
Columns 9 through I 6 

ZI 

F 8.3 format 
Columns 17 through 2^ 

XJ 

F^.3 format 
Columns 25 through 32 

YJ 

P ^.3 format 
Columns 33 through hO 

ZJ 

P 8.5 format 
Columns Ul through W 

XK 

F 8.5 format 
Columns k^ through 56 

YK 

F 8.3 format 
Columns 57 through 6k 

ZK 

F 8.3 format 
Columns 65 thnmgh 72 
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Component shapes in the existing programs are 


Shape 1 
Shape 2 
Shape 3 
Shape 4 
Shape 5 
Shape 6 
Shape 7 
Shape 8 
Shape 9 
Shape 10 
Shape 11 
Shape 12 
Shape 13 


Discrete Mass 
Cylinder 
Truncated Cone 
Torus 

Beam (Arbitrary Cross Section) 

Circular Road (Arbitrary Cross Section) 

Spherical Segment 

Sphere 

Hemisphere 

Parallelepiped 

Swept Trapezoidal Panel 

Symmetric Swept Trapezoidal Panels 

Curved Rectangular Panel 


This program can easily be modified to include additional shapes. 

Frequently where precise weights of components are knn«r, ^ 
convenient to input this weieht rather than T are known it is more 

to be calculated. The value of Thll he have 

inputing Rho as density Qb/cfi;? a T limiting value for 

weight Dbs). I„ eomfoLL “efekeArar'Sir;.'ll“‘?h-'' " 

the t,„po„ent »at be lnh„t. Thn'etln? n^«hT»^ the total weight el 

shapes can be determined fro. the "shape data Input InTtr^ItljM"!"* 

with h"!l»:*oTvo?3s'' mfconpon'e'ri‘"%“'' T* '’topertles of a component 
hollows or voldrjnaude™ as s'urLjL sT' 
input as standard shapes having negative valuer fnr^ 

program will compute the actua! Sffi^eInie"of gr^^"^^ or density. The 
products of inertia of the component. ® ^ moments and 


the «:\'en'd'f cX^^SrK:ra^:”p-'‘';e'd1es^:c^|L^^ 

can be found by "^fS^rjo tL ^Lr^ °"“ted but 

ua Dy lererring to the computer program in Appendix B. 

Note that component axes have been located so the^ v -t ^ 

axis and the y and z axes either coincide with, or Ire paSLj'to orl 

-r:-t,Tra“ti:^pa?± 
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SHAPE DATA INPUT INSTRUCTION 


SHAPE 1 
DISCRETE MASS 


Mass 



Center of gravity 


IXXCG = A 
lYYCG = B 
IZZCG = C 


^ X 


Input Data: Item, Description, Shape. RHO (wt,) 
A(lxx), B(lyy), C(lzz) 

XI, YI, ZI, XJ, YJ, ZJ, XK, YK, ZK 


OPTIONS AVAILABLE 

No.l Input data just as indicated above. 

No. 2 If negligible, inertias A(lxx), B(Ijt) and C(lzz) may be omitted, 
in vhich case points j and k are not required. 

NOTES: Never locate point j at the system origin and if input, 2 

A(Ixx), B(lyy) and C(lzz) must be inertias about principal axes in slug -ft. 

y SHAPE 2 




Input Data: Item, Description, Shape, RHO (density) 
A, B, C (option No. 2 only) 

XI, YI, ZI, XJ, YJ, ZJ 


OPTIONS AVAILABLE 

Z‘c.1 Input data just as indicated above. 

::o.2 Cylinder may be segmented requiring a C value be input. 

;;o.3 Total weight may be input for RHO in both previous 

options but only if it is input more than .k pounds. 

Xo.b Input A and B equal and total weight for RHO 

and program treats shape as a thin-wall cylinder. 

POLES: Point k is not required with any option and density 

as such mush never be input more than .U pounds per cu. in. 


POOR QUAUiX 



SHAPE DATA INPUT INSTRUCTION 



Input Data: Item, Description, Shape, RHO (density ), A, B,C^D 
XI, YI, ZI, XJ, YJ, ZJ 

OPTIONS AVAILABLE 

No.l Input data just as indicated above 
No. 2 A total weight may be input for RHO but only if 
it is more than ,h pounds, 

NOTES: The values of A minus C and or B minus D must never equal zero. 

Point i is always at the cones Larger end and point 
k is nat required . 

Density as such must never be input greater than .2; Ibs./cu.in. 



XI, YI, ZI, XJ, yj, ZJ 

NOTES: In this case program assiimes all muss 

is concentrated midway between inner und outer surfaces. 

Point i is always at the cone's la7*,‘‘rr end and point k is not re- 
qtiilred . 



SHAPE DATA INPUT INSTRUCTION 


SHAPE h 
TORUS 


y 




X 


IXXCG = Y11-Y12 
TYYCG = X11-X12 
IZZCG = lYYCG 


Input Data: Item, Description, Shape, RHO (density), A, D 
XI, YI, ZI, XJ, YJ, ZJ 


NOTE: Point is not required. 


SHAPE 5 



Input Data: Item, Description, Shape, RHO (density), A (area), 

3(lyy), C(lzz) (Area moment of inertia in inches^) 
XI, YI, ZI, XJ, YJ, ZJ, XK, YK, ZK 


NOTES: If beam is a body of revolution about the x (centroidal) axis, point k 

is not required. lyy and Izz are area moments of inertia about 
principal axes of the beam cross section taken in a plane normal to 
X axis • 
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SHAPE DATA INPUT INSTRUCTION 


SHAPE 6 
CIRCULAR ROD 



NOTES: 

Point "k" is not required and weight Instead of density is the required 
input for RHO. 

The solution for this shape is approximate in that it may be as much as 
one $ less than correct if rod dimensions t and w are as much as 2% of 
dimension A. (Error incurred tends to increase as this ^ Increases) 



No.l Input data just as indicated above 

No. 2 Total weight may be input for RKO if more than ,k lb. 


NOTES : 

Point k is not required. When computing as a solid spherical segment, S dim. 
will become equal the distance between points "i" and "j". Density as such 
must never be input greater than .4 Ibs./cu. in. 
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SHAPE DATA INPUT INSTRUCTION 


y 



SHAPE 8 
SPHERE 


IXXCG = XI1-XI2 
lYYCG = IXXCG 
IZZCG = IXXCG 


Input Data: Item, Description, Shape, RHO (density), A,B 

XI, YI, ZI 


NOTE: Points ”j" and are not required. If B = A the program selects 

thin-wall equations; therefore, weight instead of density should be 
input for RHO. 


y 



XI, YI, ZI, XJ, YJ, ZJ 

!;CTES: Point is not required. If A = Lgth (i,j) the program selects thin- 

wall equations; therefore, weight insted of density should be input for 
RHO. 


PA-GB ^ 
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SHAPE DATA INPUT INSTRUCTION 


SHAPE 10 
PARALLELEPIPED 



Point k 

(Can lie anywhere 
in the xy plane) 


,(^- 1 + 

lyi-jj 

r I J I ,'T 


A 

F 

p 



Inside Void 


4 

k --- c 

(Must be centrally located 
inside parallelpiped ) 


IXXCG = XI1-XI2 
lYYCG = XI3-XI4 
IZZCG = XI5-XI6 


Item, Description, Shape, 
RHO (density), A, B,C,D,F 

XI,YI,ZI,XJ,YJ,ZJ,XK,YK,ZK 


OPTIONS AVAILABLE 

No.l Input data just as indicated above. 

No. 2 A total weight may be input for RHO 

but only if it is more than .L pounds. 

No . 3 Input D equal A and the program 

selects thin-wall equations therefore 
total weight must be input for RHO. 


NOTE: 


Density as such must never be input greater than .L Ibs./cu.in. 
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SHAPE DATA INPUT INSTRUCTION 
SHAPE 11 

SWEPT TRAPEZOIDAL PANEL 
(THICK WALL) 

Input Data; 

Item, Description, Shape, RHO (density), A, B,±C,F 

XI,YI,ZI,XJ,yj,ZJ,XK,YK,ZK 

(Assign i to C according to sweep of panel) 


Formulas designed to 
facilitate input data 
determination : 


= £(Bl 


= 

3( 


B-F 

B+ZF) 
B+F) 


^ ^ L^th(B+£F) 


3 (B+F) 

OPTIONS AFAILABLE 
j- No.l Input data just as 
NOTES : ' indicated above 

Dimension F must never be more than 9??^ of B. No. 2 Total weight may be 
- ensity iuch r.io.st ncv^i. input greater input for RHO if 

than .k Ib./cu.in. more than lbs. 

SHAPE 11 (SPECIAL CASE) 

UNSWEPT TRAPEZOIDAL PANEL 
(THICK WALL) 

Input Data: 

Item, Description, Shape, 
RHO (density) , A, B,F 


XI,YI,ZI,XJ,YJ,ZJ,XK,YK,ZK 



Point k (Can lie any- 
where on the y axis 
excepl at point i ) 


OPTIONS AVAILABLE 
No.l Input data just as 
indicated above. 

No. 2 Total weight may be 
input for RHO if 
more than .4 lbs. 

NOTES : 

Dimension F must never be more than Q8^ of E. 
Density as such must never he input greater 
than .4 Ibs./cu.in. 21 







SHAPE DATA INPUT INSTRUCTION 
SHAPE 12 

SYMMETRIC SWEPT TRAPEZOIDAL PANELS (THICK WALL) 


P6int k 

(Can lie anyvhere 
in the xy plane) 


Input Data: 

Item, Description, Shape, RHO (density ), A, Bi--C,D,F 

XI,YI,ZI.XJ,yj,ZJ,XK,YK,ZK 

(Assign • to C according to sweep of panel) 

— - Axis of Symmetry 


r I u 




EXAMPLE PROBLEMS 


IVo example problems are presented In order to show required input data. 
These problems were selected due to the simple calculations involved and thus 
could be checked by hand calculations . 


Example Problem 1 


Problem 1 was taken from reference 3. It is a cylinder skewed in the 
y -2 plane. The moments of inertia about the system origins are given in this 
reference and are used for comparison in this paper. It should be noted that 
"g” (acceleration of gravity) was taken to be 32.0 ft/sec^ in this reference 
rather than 32.2 ft/sec (386 in./sec^). Depending on the value of ”g” selected 
by the program user, the term 


**cons** has to be changed accordingly. 


The axis through the center of the component must always be the x-axis 
for the computer program. Therefore, the y 3 axis of this problem corresponds 
with the x-axis of the computer program; the yi axis corresponds with the y-axis 
and the y2 axis corresponds with the z-axis . 

Find moments of inertia I*ii> 1*22, 1*33 which corresponds to Ixxo> 
lyyQf Izzo coordinate system. 

Given: W * 1 slug = 32 lb 

R = 24 inches 
= 36 inches 

e = cos 3/5 = 53° 8' 



Solution ; (a) Direction coelnes for 
the transformation are 


‘IJ 


1 

5 


0 - ^ 


0 

k 

5 

2 

5-J 


And ntsoerlced values of the 

' p» O 0 

hj “ 


■■ij 


are 


Component 
Moment of 
Inertia 


ORIGINAL PAGE IS 
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0 
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The results of problem 1 as computed in reference 5 are now given. Hence, 
from the transformation equations 

Gy stem Moments of Inertia 

^ 12 ^ 12^22 ^ 13 ^ 15^55 ^^ 12 ®' 13^25 ^^ n ®’ 12-12 

I^x = U + 0+ 0 + 0 + 0 + 0 = “= sl’og (Ixx = 

-22^ ^21®’21^11 ^22^22^22 ^25°'23*53 ^^22°’25^25 ^^21^25^15 ^®'21^22‘12 




° * (f)'' ^ ‘(sf ° ° ° ’ i “ * ‘2^' 


Gi " ^51^31^11 ^52^32^22 ^ ^33^55^33 ^ ^^32^33^25 ^^31^33^13 " ^^3l''32^12 


-ZZ “ ^ 


h 2 52 

" -5 '*5 


r^O ^ p 

+ G-rO=^ = 3.20 sl'Og ft 


?ni5 problem is now computed using Computer Program roCTC2. Sho"*^. below are 
:hc coordinates of points i, j, and ic for problem 1. 

V > ^ 

< --: = 2 

v^z. = C 


/ iZ = 3'C sir. 53“':® ' =25.; 

> 

i^Z . = 3o cos 55"^® ' = 21.; 


Since not recuirea lor cv^i nc 


.-lerence us now maae to tne necessar;/ data cards to compute this problem. 
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ITEM NO. 


lii 


DESCRIPTION 


IHDER-TEST 


WT. OR 
DENSITY 



















We now have the results of problem 1 using Conrputer Program 

(Not* that th« gravlUtloMi consunt ustd hert was 32 . 166 ** ratiker than 32.0 ustd in rtf.) 


ITIPI 

I 


ottCBitrio* 

■I »i 

••»§« ••••• 


iMsn •*•0 A a c 

I .i20»9Of»Cttl M.690 «•#•! 9.999 

II rj IJ i« 

9t999 9.991 I9.a99 tl.a«9 9«999 


Q 
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Example Problem 2 

Problem 2 was taken from reference 2. In this space structure, 
the weight has oeen ’lumped’ or concentrated at the joints. In the 
program being presented, this is not required but is used here only 
for illustration. To compare moment of inertia, the numbers given 
in reference 2 should be converted to slug-ft^. 


IXX 


22038.464 

4.608E3 


4.783 slug-ft^ 


IZZ 


9000.00 

4.0608F3 


1.953 slug-ft2 
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ExompU Problem 2 





INPUT 

MEMBER AREAS -OOlm* 

WEIGHT AT EACH JOINT - (5 Q lb 
MODULUS OF Elasticity - lo^pt, 

SEE the first two PAGES OF THE SANV>LE 
problem FOR details OF THE INPUT 


JPiMr 

- JOtkf II 
1 0. 

t to.oouo« 

y lo.opooc 

* 0 . 

« 0 . 

ft 0. 

V io.»oooe 

I 10. DOOM 

ft 10.00000 

10 O. 

11 0 . 

11 lo.ooooo 

l> f.OOOOO 


u 

0 . 

0 . 

10.00000 

to.ooeoo 

10.00009 

0. 

0. 

10.00000 

10.00000 

10.00000 

0 . 

0 . 

9.09000 


II 

0 . 

0 . 

0 . 

0 . 

10.00000 

10.00000 

10.00000 

10.00000 

20.00000 

20.00000 

20.00000 

20.00000 

10.00000 


The results of problem 2 as computed in reference 6 are now given. 


fi lAL 105.000 iftl.OOO 105.000 0. I. ^ 

..UtlftiUlf-Ol li.ht I* 5.000. V. 5.000* 1* 11.550 _ . . . - - 

HONCNT'. ur lOfOTU ftOOV* CINttO Of tfOISMT k 

111 * 220 i« ftftJ. ivft. 22 OI 0 .fft 2 . Ill* 0000.000* lift* >0. i |l|* O.OOOi IVI* 0.000 

;;OTE: Inertias are given in pound inches sr,uar<.-G ace be converted zc compare with results progra=. yiel-a. 


This problem is now computed using Computer Program ND0702 
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Item No.l 
(1.5 It. total 
wt. ) 


Sta. 35.0 


EXAMPLE PROBLEM 3 
A COMPOSITE STRUCTURE 

( ^ SCALE ) 


^ Item No. 2 

/ ( 1.5 in. thick) 


Sta. 25.0 


Sta. 0.0 

Item IIo. k- 
(o- thick) 


Sta. 



Item No. 6 ^ 

(Cutout in Cylinder) II 


Sta. 11.0 


Item No. 


NOTES: 

Points i and j shown on the struc- 
ture Eire those assumed when data 
was input to the inertia program. 
Also, some salient dimensions are 
included on the components to 
assist in relating component data 
to the computer input. Except 
for i’tet^^s -o* 1 end No. 7 de^-sit;' 
.'HDut for RHO. 


Item No. 7 

(5 Ih. wiring etc. evenly 
distributed inside cylinder) 


Item No. 5 


JL. 


0 in. 


The computer output relative to example problem 3 follows and 

the , "*he *^oinp^ter calculated component data, 

the soTiined data and lastly, the inertias about the system principal 
axes and the individual axis locations. 
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CONVERSION TO INTERNATIONAL SYSTEM OF UNITS 


•NVERT FROM 

TO 

MULTIPLY BY 

INCHES 
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.000 6^5 i 6 o 

FEET 
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,'^Oh 800 
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PObTIDS / inches^ 

]Q LOG RAI VME TE R^ 

27 679-905 • 
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kilogram/neter^ 
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SLUC- 

KILO'G?Ai''iS 

Ik. 595 9029 

2 
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2 
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CONCLUSIONS 

A computer program for determining mass properties of a 
rigid structure is presented. The structure is broken down 
into preselected shapes with known properties, and input data 
are supplied to completely describe each shape. For compli- 
cated structures skewed in space, this program offers a practical 
solution to a tedious and time— consuming task. It is also prac- 
tical to use this program for problems that involve repetitious 
cr lengthy calculations. 
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APPENDIX A 


DIRECTION COSINE DERIVATION 

The direction cosines are used to transform the properties 
from the component coordinate system x,y, z to system coordinates 
X.X, Z. Three points (i, ,i, k) shown in the figure below, define 

the vectors V_ and V ^ 

1 2 * 



Vector is arbitrarily/ selected to be coincident with tne 
X axis. A unit vector on tnis axis can be written 


1 


' V 


'l! 

Taking the vector cross product of with ’.V and dividing 
one resulting magnituie gives a unit vector on the z axis 


triiy, a uni* 


on the 


axis IS fcuni 


] X 


The direction cosines are tne X, Z, ccriponents 
vectors cn the x. y. z axes. 


The direction cosine for the x axis Is written as 


LX = (XJ 

-XI 

)/LGTH 

MX * (YJ 

-YI 

)/LGTH 

NX = (ZJ 

-ZI 

)/LGTH 


where the length Is 

LGTH * ((XJ -XI +(YJ -YT +(ZJ -ZI 


The vector V 2 can be written 

T1 = XK -XI 
T2 = YK -YI 
T3 = ZK -ZI 


A vector on the z axis is found by taking the vector cross product 
of \ and V 2 

LZ = MX*T3-T2*NX 
MZ = NX*T1-T3*LX 
NZ * T2*LX-T1*MX 

The length Is 

T4 = (L2^ + MZ^ + NZ^)** 

Normalizing to get a unit vector 

LZ-LZ/T4; HZ-^Z/T4; NZ-^NZ/T4 

The unit vector on the y axis has a^magnitude of one and Is determined 
by the vector cross product and 

LY = MZ*NX-MX*NZ. 

MY = NZ*LX-NX*LZ 
NY = MX*LZ-U*MZ 


Writing the nine terms In matrix form, we get 


fOC]= 


LX 

MX 

NX 

LY 

MY 

NY 

LZ 

MZ 

NZ 
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APPENDIX 6 


PROG«_AM name ( IMPUT .UUTPIIT t T A PP 5 = I MPUT , T APF 6=01)1 PUT ) 

C LATEST DIRCQS- IN MAIN DECK AMD 13 SHAPES 12 / 70/12 


- PJMENJ5L0!i__II£M( 2 i )0 ) , SHAPE ( ?oo ) , Rmo ( ?oo ) , a ( 20 ( 1 ) , R ( ?nn ) , c ( ?nn ) , 

1 D^^ 00 ),F(^ 00 )^ K I ( 200 ) , Y I ( ?nO ) ,z I ( ROO ) , XJ ( 200 ) T YJ ( 200 ) . z ) ( ? 

POO ) ,XK ( 2 OOL 1 YK ( 200 ) ,ZK ( 200 ) , I XXCO ( 200 ) , I YYCG ( POO ) , T ZZCG ( 300 ) , T X Yf.r, 
3( POO) »1YZCG(200) tI_XZCG( 200) .XCG( 20 0) ,YCG(200) ,ZCG( 200) »XL( POO) ,Y| ( 
A200) tZL ( 200) , I XXCni 200 ) , I YYCP ( 200 ) , 1 Z ZCn ( POO ) , T XYCO ( poO ) . I XZC'K PO-^ 
S ) , IYZC 0 ( 200 ) ,DES( 3, POO ) ,w ( 200 ) .n( 200 ) , p ( 200 ) , ARR ( 3 , 0 ) , p ( 3 ) . r ea / ^ ) 
real IKX.lYY^,IZZ»IXXCG,IYYCO,IZZrr-,IXYCG,IXZCr,, 

1 IVZCG»LX ,*iXiiNX « LY *MY .MY ,LZ .mz ,m7 , TXY, IXZ» TYZ.LfiTH. TXXCn, IvyC'’, 

PIZZCn.IXYCn,IXZCn,IYZCD.IXXn,IYYn,i77fi,KI,KII,Kin.IXYK.IYYrOP 

IMTEGER SHAPF 
1=0 

1010 1 = 1+1 
106 FnP>-AT(I5)_ 

R E A 0 ( 5 , 10 l) ITF:V1(I),nPS(l,I),.0Ps(7,T),SHAPP(I).QMn(I),A(il.M(M. 
lC(I).0(I)^jF(I). XI(I).YI(I),7I(T),yJ(J),yl(l)^7l(J)^y^^(|l 

PYK ( r ) ,Z*< { I ) 

10 1 FORMAT n 3 ,2A9 , I 2,F 9 ,SP ° /QP ) 

I F ( PPF ( 6 ) ) 1000 . ] O 1 0 
loro IMAX=I -1 

OR I TP (M. 271 ) 

P71 FORMAT)] MO* IMP! IT data LIST-'^ p = ).■''•>*// ) 

WR! TP ( 103 ) ( I TPM ( I ) . ( 1 . r ) . ( 7, T ) . vhAPP ( I ) , ( T ) , A ( T ) , - f I > . 

1C(I).0(I),F(I), XI(I),vj(j),7j([)^y|(j)^Y.)(T),7J(]t.xk(T'. 

PYK ( I ) .ZK ( I ) , I = 1 . f.,AX ) 

103 FOR '-'AT (2XX--ITP'-. OP I -^T I o- 

1 f 

P -■/!‘'.''X,2AR.IS,lX.Pl^,c,‘,,-1-..:i/;: 

3 zi x,l YJ 

A YK 

1- R ! TP ( ^ , 261 ) 

PAl P(lR A T ( 3X9MC QP'PDimPm r PX4MI I at A P y Ami IST‘^ 

HR I TF ( 6 . 250 ) 

PSD POS p A T ( 1 '-'0=:= ITPM OPSCPlMTir.. 

1 Imyco IZZCO XC'^ro vror.- 

p ixzr,,.::^// ) 

T-: = 0. 

X n V = n , 

Yi-',nw=o. 

7 MOM='''. 


S A P (“ i? '1 

'<1 V T 

/.I 

' V xS^<=iC|_riGi// ) 

'-T Tvvr. 

7 C r. c p T y V c ‘ ' t y 7 f- 
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.. . ... 



m«o. . - — - - - 

Dlyjlo, 

IXZ«0, 

IYZ«0, 

ij(xo»o. - 

IYV0»0, 

IZZO»0. 

PI»3.iai59S 

COMi«4632, 

D050 I«1,IHAX 
XLn)»0. 

YLn)«0, 

ZL(I)»0, 

IXYCG(I)«0, 

I*ZCG(I)«0, 

IY2CG(I)«0, 

LGTH«so«?T( fxj?i)-xi(n)**2t£Yjm»vitn)t*2*(zjin»zi(in<»*2) 
ITt»<P«5H*PEtn 

GO TO (1 ,2,3, Cl, 5,6,7. !0»n*12»l3),ITEMp 
c OISCI^ETE '’ASS *SHAPE I* 

C point J '’JST be selected ANVfcHEPE On Th^ x AXIS FOP DIP, COSINES 

1 --(nBPHon) 

i*xcG(i)»Acn*coNS 
IvYCGdlaBt n*C0NS 
TZZCG(I)BC(n*CO'.'5 

o,).o‘».(APS(Tjnn,Nt,o,),op, fABsczjcni.NE.o, 

1)) G-1 TO 7jil 

UGThsI 1 ), 

*J(T)s1,l*'Cin)*10, 

IF(Xj(I}.F'i,*I(n)xJ(l)*Xl(I) + lO, 

7 Oil 

Ml ir. fru 


c CtlI^'OEw *S'*APE 2* 

2 IPfE>(n.E3,Afn5GO TP 53 

TFfprn.LT.o.i'^fns^'. 

I^CC(I).LT,C'.ICm«n. 

v^'LCxi. sPT*f Afn**2-Rf n**?i*rLGTH.cmi 
*f 1 )«PHPf I )*i/OLCVL 
IF fp«nn).5T.o, tn KfnsPHorn 
pup f Dsn f T I /vOlCVU 

i*xcorn«f, 5 *-‘M)*fA(n**?*flfn** 2 n 

ivvcGfn*.?5*'^m*rA(n**2*Bfn**2)A<«ifn*fLGTH**3,rfn**3> 
*/f 12,*flGT-<*Cf Tin 
izzcGrn«ivvccf n 

KLfllBLGTH/P, 

GO T r. 0 


original PA^ 

OF POOR 
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cylinder (THIN NALL) *SHAPE Ia 


33 '*(I)«RHD(n 
XL(I)*LGTH/2, 
ixxcG(n«(»MO(n*Af n**2) 

lY''CGn)«,5*RW0(l)*AfI)**2^RM0n)*(LGTH**3»C(I)**3J/ 

2((LGTH.cn))*12,01 

IZZCG(n»lvyCG(I) 

GH TO ftO 


TRJMC4TE0 CQNt *SHAPE 3* (REF--R* HULL) 

3 IP(B{I),EQ.A(I))G0 to 30 
VOL ■ 

Sl.O«72«LGTw*(A(n**2fACI)*C(n+CCI)**2«0(I)»*2«B(I)*O(n*D(I)**2) 

• fT) ■ PHon)*vOL 
iPfRHonj.GT.o.tt) *(n«H-nfn 
RHO(I) B »>tn/vOL 

*L(I)« (LGTH**2J*(An)**2 + ?,*4(n*cn)*3,*C(I)**2»«rn**?*2,*B(I3* 
TD(I)-3,*O(I)**2)#,26ie/V0L 

<I s.95N3.l«l593*R^^nfI)*( (C (l)**B-A(n**5l*LGTM/(Cf 1)*4(I) 1. 

$ (0( T ) rf)t*5)*LGTH/(n(n,0(i))) 

•'II*3.iai6*RH0(I)*LGTM**3*( f A£i)*,2-B(n-»*21/1,0-,5N( (4(13* 
TC(I))*A(7)»(B(I)»D(I)3*B(n)t,2*(f4Cn»C(I33**2*(B(I3«0(n3**2>3 
<111*1,0 u72*(LGTm3*(A(I)**2t*(I)*CCI3tC(I)**2»B(I3**2-B(I3*0(I)» 
?on3**23*RNQ(n*xLtn*)»2 
TxxCGn)* 2 , 0 *(i<T ) 

! ' ■ C ■. ( i 3 s< 7 ♦*. I I-K I 1 1 

IZZCGf n*I VYCGt 1 3 

GO TO ^0 


TijjvCATE^ CP‘Jfc (ThI» i^ALLl *So*PE j# 

33 VOLS^fI3*S3RTff*(I3«C(733**P*LGTN**23*3,lU159*fA(T3*C(I33 
-■fi) s vot*R^Ofn 

*LfT3sLGTH/3,*tf?,#cfT)*Am3/rrni*A(nn 

iVYCr.f nsfOHOf A(l r n**23*3Hnn3*L5T-**?/l<».*M ., 2 ..Af 

iTi*cfn/(Afr3>cn)3**2n . t v 

TZZCOf I)«I vyCG(T) 

I**CGn3 = (‘?Hnm/2,*fAfn**2*CfM**233 

GO TO of) 


TOP'.'S *S-'APF <j* 

y A ( T ) BA f I 3 .LGTh 

fOm.LT.O.lOf 13*0, 

'''^-1*2,*PI a*2*LGTw** 2*4 ( j 3 

'''^L?*2,*PI a*2*0( I 3 **2*a n 1 

AC v^L *VOL 1 •VOL? 

A ( T ) sfiMfi f T )* AC VOL 
X“lsfPw^fT3*V0l13 
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Xl !■, 125«X‘<1 * I ) Jt2*5,*tCTR**2J 

Xy2«(».H0XI)*VQ12)-- 

XI2«.l25*X«2*(<l.*AtI)**Pf5,*f>(n**2) ' 

IvvCGn)«fXIl-xi2) 

I72CGtI)slvyCGtn 

vTl«,25*X'^!*f«,*A(n**2f3,*LCTH**23 
YI?*,25*XM2 *(<j,*a n ) 

lxxcGm*Yll-YI2 

GO TO bO 

B£*« (ARBITPARy CROSS SFCTIOM) *SHAP£ S* 

5 XLtn«LGTH/2, 

ixxcGn)*(RHO(n*LSTH*{Bfi)fC(n)) 

IZZCGfI)sfPM0(I)*(C(I)*LGTM + ,0833#An)*L5TW*)»3)) 

YDL*A ( I ) *LGTH 

f 1 3 *rhd ( n * vhl 

GD Tn tfO 

CIRCULAR POD (APBITPARv CROSS StC^lOM *SHAPfc 6* 

6 *n )sPMon ) 

RSCG=fAfn*A2*c(n*Bnn/fAfr)ACfn) 

1 XXCG f I ) sRSCG**?*k. ( 1 3 
TVVCGfI3 = ,5*RSCr,**?**. fl3 
T/ZCGf T 3 si Y vCG t n 
-■jr t: ,0 

SPHtPTCAL S^G“FUT *ShAPE 7* 

7 P f 1 3sc n 3*Bn 3 
''fl)=r fI3-LGT>- 
G SLGT«-.a(I3 

VOL1sl,r>o7?*Lr,TH**?*f5.*rfT3«LGT^3 
y''U?= 1 . ''‘J7?*G **?*f3,*Ff n-G 3 

ACw"'Lsv'^L 1 •vOL2 
• f T 3 s AT *PHO ( T 3 
iPf5'-.'3n3,GT,o,P3 wnsPHOfn 
( 13 S • n 3 /AC VDL 

*!?AP3«,7SA(2,*cn3»LGT^-3*A2/f3,*C(T3-LGT-«3 
XWAP?s,7S*(?,*F (13-r, 3**2/f3,*PfT3-G 3 

XI y s ( y-iA P 1 * vDL 1 -*B AR? A vOl,?1 /AC vOl 
»Lm«XLT -nfi3 
X« 1 r vOL 1 *B>^n { 7 ) 

X ? s V ! ? A 9 M n n 3 

*naf?,ALGy"i*XMl/(3,AC (I 3-LGTM3 lAfCf T 3 AA?,,75*CtI3*LGTMA,15*LGT 

1 AA?3 


OMGINALPA^ 
OF POOR QUAUTX 



XI2«(2,aG *X.M2/t3,*Ff n-G ) 5 1 { F U ) dft2«,75*r ( I) *C ♦,15*C 

1 **?1 

lXXCG(n*(XI!«Xl2) 

TEMPi«,05256*HH0(n*n5,*c m**<»*c(i)»io,*c(n**?*C(n**3*5.*c(n 
U*5) + .209<*t|*«HQ(I)*t5,*C(n**2*C{n#*3-3,*C(IJ**5) 
TEMP2»,05236*RHO(I)*t 15,*C (!)*»« oDf I). 10,*C(n**2*Orn**3*’l,*D(n 
l**53 + ,209«a*RH0{n*tS,*C ( I3«(>2*Dn 3 #i»5«3^*0 f IJ **51 
TE>'P3s,0S236*RH0(I)*M5.*F f 1 1 •« «*F ( II • 1 0 , *F ( 1 ) ** 2 *F f 11 **3*3 . *F (1 1 
l»*5) + .209««*RHO(I)*t5,*F tn**2*Fm**3»3,*F(!)**51 

TE''P9«,05236*RND (lD*n5.*Fn)**0*0(I)«10,*Pn)**2*OrTl**3*i,*0(I) 
1 **51 ♦,20 90U#Rwon l*f5,*Ptn**2*n(l)**3,3,*!)fi)**si 
*CTEmp«( iTEMPl.TE»*P2)*(Tf MPi.TF^-PUl 3-(TH*XLT **21 
IvrccniBACTEMP 
IZZCGdUlYYCGm 
GO TO Ml 

SPhERF *S.hAPE 8* 

9 LGTM»4fi) 

XJ(I Isxi ( I }*LGT h 

YjfllpvTm 

zjiii»zim 

lF(?f n.EO.AmiGn TO lu 
V0t1=u,l98T91*An)**3 
188791 *B f 11**3 
AC v'TLrvnLi .VTL 2 
- f T lsR«n n 1 *AC vOl 
t « 1 s f ur ( T 1 * vriL 1 1 

X“2sfw-n(ii*ynL?l 
* 1 1 5 ( , u*X‘' 1 *LGTh* *2 1 
xT2*f ,4i*x‘^2*Rm**21 
lXXCGril3fxTl«»I21 
TvycGfTlsTxxCGfll 
I Z/CG f I 1 sT XXCG f I 1 
C-n y' e*r 

fT»4lv *ALll *SH4Pr a* 

3u 

lx*Cr;fTlsf,t)h7*>J-*^ril«4rn.*?) 

T YVTG f T isl * *CG f ! 1 
T7ZC-f naTxxCGf! 3 
X L ( T ) - ^ r 

%’'■ T "I *>n 

-If '‘T SP»~F «F *S-'APF 9* 

9 I ^ ( LGT^-,f 5, A f T 1 ifin T? 3S 
TffAdl.LT.O.lAfnan, 

'•''■■L’*2.^9i*39*i,GT*~**3 

y '-';.2»2,f'’9‘i39* A (11**3 
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ACVOLUVQLj-VDLZ 
•<(npRHr)f n*ACVDL 
xPAR1»,375#LGTH 
XBAP2«,175*4 (I) 

XLniifX6ARl*VDL.l-XBAR2*V0L2)/ACv0L 
Xmi«(RhD(I J*VDU) 

XM2,CRHDtn*VQl,2J 

xnsC.<i*X*^l*LGTH«*2) 

XJj2iX.«*X'^?*An)*<»2) 

IXXCCt J)a(XIl.XI2J 
XI 3 b ( ,26*XMi *lGTM**2) 

Xiaa( .26 «Xm 2*A(I )**2) 
IYVCCtnifX13*Xl<J) 

I7ZCGtn»iyycC(l) 

GO TO 6ft 

HE“ISPHFRE i^ALL) iSHAPE 9* 

35 ■^msRHCnj 
XL f n*LCTM/?, 

IXXCG(ns,666*RHn{I)ALGTH**2 
lYyCGm«.«166*Rwri(n*LGTH**2 
IZZCGtDslYYCGfn 
GO T<i 6 0 

PARALLFLEPTPEn *S-APr 1^* 

10 TP f o, ( I 1 ,t3, A n n GO TO 36 
IFtOtD.LT.ft.lOtnsO, 

VOLI =LGT*’*a f n *A ( I) 
v'^LPsc n ) *p t n *0 fi ) 

AC^OL ryni_ i ■V.0L2 
‘fTT=‘«-' fT)*Ar 

j P ( a M p f T 1 . T , 0 , (j 1 . n ) s a « n r T I 

p-'^d) s »r!J/ACvnt 
XL ( I )SL r.T-/?. 

XM] syr>L 1 *Buri( I ) 

*»i’ 2 synL f I ) 

»nsf.f'6 5333*x“l*fRrn**2yAfl)**2)T 

xi2Br,OB^333*x«2*(Pfi)**2^0fT)**2)^ 
Tmrcr, (Dsfyn-xT?! 
*T3sf,oft3535*yMj*(L(-.Tw**2t4n)**2n 
>Tasf,0«33?3*X''?*fCm**?*Om**?n 
ivvcGfnsrxTi-xiu) 

*TSsf,o«3533*X''l*fLGTH**? + a(T)**2)i 
xT6sf,.'TB3333*x“2*fC(n**2*r(n**2n 

TZZCGf I )s(y T5-XT6) 

G!'; T-' 6r 






0^ 



P*B*LLELEI>IPKD fTHjN wiLLJ *8H*Pt 10# 

5fe xL(I)*LCTm/2, 

TE‘^Pl«(LGTH*B[n*A(in 

TE'-PJ^fLGT^^^BfD+Bd)** m#LGTH*Af n J 

IXXCGn)!iC,083133*PH0{n* fB(n**2 + A(n#*2)#(RH0n)/6,)*(TEHPl# 

1 (B(io#*(ii)/Te«p2)) 

iyrcGfn«( ,oe3333*PHO(n*fLGTH**2#A(n**2itfRHocn/6,)*nc»*pi* 

1 (L6Tw + 4m)/TEMP2J) 

IZZCG(n*C,0B3333*PH0(n*CLGTN**2*B( n**2U(RH0n)/6,)*(TE‘‘Pl* 

1 tLGTH + Bn n/TF**B2n 
GO TO 60 

5«tPT trapezoidal PA»JEU (THICK oALL) tSKAPE ll* (REP«-«R. -•ULL) 

11 ••ti)*ArnirRMon)*(LGTH*(H(nfF (in/2.) 
iP(OH0f n.GT.o.u) •.(nsR^ntn 
RH 0 (I)s(-.n)) /(((Bcn*^ (1) 1/2,0 l*LGTHiA( in 
xL(i)«LGTMA(B(n+2,o*F cin/(5.o*(B(mFn)n 
XT S Fn)*LGTH/(B(I)*FCin 
T^TAIJ * Cdl/LGTH 

FETA^ * fF(I)/2,- 9dl/2,f Cdll/LGTM 
AFTA^ * (Cm-Fd5/2,0*Bm/2,0)/LGTH 
thPB«lGTh*xT 

BCG*C(KFTAV4.AFTAV)/2,l*(Twt?B.)(i.d)l 
xBA^siF’T Av.trr AK. 

IVVCGdl* H-i,id)*{(XPA'.«fn**3*(TMHB**2-XTA*2n/2<l.* XOA^AAdl* 

2 (THPfi**(j.wT#*u)/u, )•( tTHf<g,xLdn*«2)**'d) 

XFRXX* xPAMAAfl)* ( P C G * A 2 * ( T HP B * ♦ 2- X T * *2 J / 2 , -BCG* ( FE T A v + A F T A ‘v' 1 # 
l(THQO**j-*T**3)/3,♦(ftTA^-♦APTA^3**2A(THR«**^|•XT**a)/16,) 

IxxCGdls RHOdl* txPA'-**3 *Adl*(THSH**<j«xT**a)/aB, AXPASAAfn**! 
I* d«P“**2-xT**2l /2u , ♦xmxxi 

IZZrGd)»Ad)*RMnd;*(XPAK**3*(Ti^Bg**U-XT**a)/<i«, + XPAM*f((ThPR-XL( 

3inA*2*«CG**2)*(TH5P**p.xT**2i/2,-(2.*(T«Rg-XLdmfFETAVAAFTAVj)*R 

ucr. l*(T-0^**3-XT♦<t3)/J,♦fU, + (FtTA^♦4PTA^J^*•21*dHP^♦*«-XT**a)/l^,n 
IF ftFSfCfTll.LT, ,00n GD T C 60 
4Pxcr,sHd)/2,AARSCT\TAj*XLnn 

MsAPxcGASDRTfixHfn-xTfinAAP+fxxcTi-vTdn^^z+czKdi-zidn**?) 

CG<S>"-APXCG 

PRHsTmRPv^CTH 

S*'‘ts"*PPB/TMPB 

ppTxVKs.BAFd)*{2,*S^'--Fdn * (LGT-**ppb/3,+prb**2/ 12, > 
lxY<sTk.pR*#p*(H**2»("»Pdn**2)/2u,-PRT»VK 
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I*VCG» i *?8(IxyK*4(I)*«W0(n •»tl)AtC'(i»XLtJJ3*Cm/AB8(Cfn5 

P*XRC«.765u**BS(C(I))ZCn) _ 

IP(*BS(I*XCGn)-IvvCG(I)).LT, ,001) GD TD 20 

(2,0*IXYCGP)/ tlYYCGn)-IXXCG(I) ))/2,0 
!P(lxyC6(n,GT.IYYCGfl)) PAXRC f *8S (2,i>P*XRC ) *3, 1 « 1 «» ) * 

1 P*XPC/4BStPAXRC) 

20 T-PJxx ■IXXCG(n*C0S{P/kxRC)tft2*lYYCG(I)*3lN(P«XRC)*^2» 

. 12,{'*lXVCCP*SI*JfPAXRC)*C0S(PAXRC) . . 

IYVCG(I)» lXXCG(n*Sl^tPAXRC)*A2+ I V YC G M ) *COS (P AXPC ) * *2 ♦ 
2(2,*IXyCGP*SI»^(PAXRC)*C0SfPAXRC) ) 
ixxc&n )bt«pixx 

2t ShFTI«*BS(TA'w(PAXRC))*XL(1) 

XJ(I)» fXI (n*(L6TW-XL(I) )*XJ(I)*Xt(I))/LGTH 
YJ(T)* (YliT)*(LGTH-XLM ) )*YJ(I)AXL(n)/LGTH 

zjm* rzif n*(LGT^^-xun)+zjn)*xLnn/LGTH 
XH n«f XI ( n* (SHFTI+CGX) -XK (1 )*SHFT1)/C6S 
Yl f ! )«{YT tn*f SHFTI+CGK) -yx tn*SrtFTI)/CGK 
ZI (I)*(ZI f !)*(SHFTI«CGX) -ZKn)*SHFTn/C6K 

l6t-*s'^pt( (xj(n-xT (T ) )** 2 *(Yj(n-vi fi))** 2 *(zj(n-zi (in**2) 

XL(!) * LG^-^ 

50 TP hr> 

Sv“-fT«Tc TQiPtZOTDAL PAVELS (ThIC< -ALL)*ShaPF 12* (REF--R, mULL) 

12 IT,T« : nn))/2,n 

V'Lsf L C-T-* (ft ( n ( n ) )*A(!) 

* ( T is =- • f ! 1 

1 p ( n . . c . u) ... f nsw-Mf T 1 

/((^‘^m*Fnn )*lgtm#a nn 
XI Misup.f-*r^(n*2,o*r(i))/(3.o*(B(n*F(nn 

Tvycp.f!) s (Lr,TH**3)* (F (n**2*«,*F(i)*B(n*Btn**?)**m*PH0(i)/ 
2f f3K.*fw n*Bf in)*.s)*(A(T)**2/i2.*(r>(nA.xLn))**2)*''fn 

XT S F n ) * lGTh/ f a ( I 1 f I ) 5 

IF(B(T)..T,Ff XTsBr^*LGT^-/fR(n-F(Tn 
T‘'T»,. = Cfll/lGTi" 

FFTi' s fF'Tl/2,- Ct.in/LGTH 

AFT* s f c f T ( T 1 /2.0 tB t n /2.0) /LGT^^ 

TFfani... ’,FfT)l LGTh s P .^1 

i»I**SRfC-'>fj;*An)*(FFTA\**3-AFTAV**3)/3,0) 

<Tj X (ftfn*A(n*(BM:)(i))/(LGTH*xTn*(TM*u*rLGTw+xT-xLfnn**? 
•TTsAftST^m 

TF f B f 1 1 ,uT ,e ( T ) ) lGT-i s M ,0-F (n/otin*xT 

XTLGT- s ♦ AHSfXn 

Txxcr. f!1Sf<T*ffxTlGTr )**iJ«flfT)i*<O/2,0)-KlI*fOrT|GTH )»*2.XT**2) 

1 ♦ AM 1 **?*-. ( n /12, 

Tz/CGm « ixxcGrn ♦ ivycgm)-'. 

XL M 1 s 1 5 TM * !5 ( n 

Lr,TrsR’3Tffxirn«xinn**2TrYjm»vinn**2*(ZJfT)-ZTnn**2i 

T ' *- 
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CURVED THIM wall panel *SHAPE IJ* - 

U «(!)• 2,*L6T«*A{n*»n ) *C(I) *RHO(I) 

IP<RHO<I),6T,0,a) w(n«RW0(I) 

HHOU) B w(I)/(J,*LCTH*Af I)*8{n*C(I)) 

»<i«j.*cn)*A(n*Bn )*l6tm*rmo(i) 

*IlBL6TH**lyJ2, 

VLn)aLCTH/2, 

IXXCC(I)bl6TH 

s *RHO(n*A(I)**5*B(i)*(2,*C(n*(2,*SlNfC(I))**2)/C(I)) 

IVYC6(nBKI 

» •(AU)**2*(C(i)»siN(c(i))wcosrc(i)))/(2,*cn)) ♦kiij 

izzccn)Bxi 

» * { A ( 1 ) **2w(C ( I )*8IN(C{ I ) )*C08(C ( ! ) )*2,*STN(C ( I ) )**2/ 

lC(I))/(2,*C(I))*xin 
GO TO bO 

BEGIm oircos 

60 IE((AB*(X<(I)),igE, O,),OR,(A8S(VK{n),v£,0,),OR, ( ABS ( ZK ( I ) ) ,NF , 0 , 
D) GO TO 90 

x«(i)Binn)*cvj(n-vi(i)) 
v<(nBvi(n + (xj(i)-xi(i)) 

ZK(1)bZI (I) 

IX((VJfl),MF,Vl(I)),OR,t'^J(n,^fe,XI(^))GO TC! 90 
'fn ( I )«LGTN*t i V i ) 

90 LX«(*J(n.XI(I))/LGTM 
-X8(vjn)-vi n))/L6TM 
N-X«(ZJ(I)«ZI(I))/LGTw 

Tl« xxn).xi(l) 

T?« VK(I).VI(1) 

T5* Z-^TD.ZKI) 
lZb >***T5»T2*wx 
‘^Zbnx*! 1«T3*lx 
'•Z»T2*lx»T 1 *MX 

T«* Sa»T(LZ*w2^ '*1**2 **^Z**2) 
lZ«LZ/Ttj$ mZbmz/t«S ^Z«^Z/T^4 
LVb mz*mx.mx*n,z 
*‘¥*N<Z*LX» iyX8LZ 
^y«Mx*LZ•LX**^Z 
E^3 niRcas 
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rotate CD’^PONfNT moment OF IN'taTIA TO S^STE^ COORDINATES 

ix*co(i)B(ntxc5(i)<t(Lx)**2*iYYCG(i)*aY)**2Aizzcf;fi)*(LZ)**2)/cnvs 
IYYCD(IT8(!XXCG(n*(*^X)*A2*IYYCC(I)*(^Y)**2*IZZCC(I)*CYZ)**2T/C0NS 
IZZCCtn*f IXXCG(n*CNXJ A*2*lYYCG(n*(NY)A*2*IZZCC(I)*(NnAA2)/C0N8 

ixYcnm* cnxxcG(n*ax)A(#^x)*iYYC6(n*av)*(MY)*izzcczn*az)A(M 
^IXZCOfl)* f(lXXCGtn«fLX)AfNX)*lYycGfl5ACLY)*(NYT*17ZCG(n*(LZ)A(N 
^lYZCOnis ((IXXCGtnAC^XlAfNXT^lYYCGfTTAfMYTACNYTAlZZCGCDATYZTArN 

iz) n /CONS 

calculate component center OF masS COORDINATES AND wRITF OJT 

XCGfl T»XI { 1 5 ♦XL (I) aLX 

vCGflTs'^I t n YXLn ) A“X 

^HiTE(*!ioJiiTE«tn .DfS(i.n,DESf2,n,-m.ixxcom,iYYco(n,TzzcD 
!(I)»XCG(I).YCG(Ii»ZCG(IT.lXYCO(l)»IYZCO{!5*TXZCOn) 
jf)'' P09‘'AT(lS,5x,2AQ,feFll,5.iJl“»,5/) 

CALCi'LATE SYSTE^' *'E1 GhT. SECOND -0-^EnT AY ORGIn AnD C.G. 


T . s T - ♦ > ri 1 

;<*'^sTxX'!‘+IY*C'^M)Y’(I5*fZCGfI)**2 + YCGCI)**2)/CO\S 
Tvv''slvvo*!YyC''(!3+*n)AfZCGfn**2*XCGfl)#*2)/CONS 
iz70s!ZZO+IZZC0n)Y*(I)*fXCG(n*A2*YCG(I5**2T/C0‘S 

j M - y 'j r* - ♦ > t n * X c G f n 
M )*VCG^ T ^ 

2^ ^ ^ ^ ( 1) ’kJCG i 1) 

rrv-piT^: 

/} A 

73A;TsZ^ n 

-wire ^ 

Sv$Tt^ DAT^ LTSTFO BF. L-V'^ (<^TsLB5i# I^’E*?TTAS*St^A*S FT Ss 
IJAPF'.', r.G.siNS. SECONH ''OmEnTsSL JG pT SOu A»E 0 ) * / / ) 

-9ITt (6. If''") Y-. , TXXD, T VYO, I7Z0,XHAR, V3AR,ZPAR 
< '- pr)Bv/T fl 1 *3'JSySl n2H^T12XUHTXX01UXaHlYV01P*P^TZZ01hXOMX3ASiaXUWVH4R 
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thasSFIR ^AIS RRORtRUES TQ SYSTEM SUM A-Nd wR.lTE QUT 

0070 1»1,1MAX : 

OELXsxCCf n«XBAR 

OELYlYCGdJ.YBAR 

DELZ»ZCC(n-ZBAR 

IXXeIxXfIXXCO(Ii-fM(no tDELYAA2 + DELZAA2) /CONS 
IVY»lYY*IYYCOtl ) + i^ f I) A (DELX*a2 + D£LZ*a 2) /CONS 
IZZ* IZZtlZZCOfDA*' t n* tDELYA*2+0ELXA*2J/C0NS 
IXV»IXY + IXVCO(n+»'fI)*DELX*DELY/CONS 
IXZ«IXZfIXZCO[n + *(I )*PELX*DELZ/CQNS 
7 0 IMZAlYZ + IYZCOdJAi^tnAnELVADELZ/CONS 
*RITfct6,200)IXX,IYY»lZZ, !XY,1X2^ lYZ 
20 0 XOUMiTngxSHIXXlbXSHIYYlSXSHlZZlSXSHIXYlSXSHIXZlSxSMlYZ/bFlB.S/i 

CO*'RUTE INERTIAS (EIGENVALUES) ABOjT PRINCIPAL, AXES ANO EACH AXIS 
OIRECTIPN cosines (EIGENVECTORS) AnO i-RITE OUT 


print 3(J0 

i<iO FCR“AT(/IX,*I\ERTIas (EIGENVALUES) AROJT SYSTEM PRTNri=»AL AXES "d 
AXIS direction CnSIvfS (EIGENVECTORS) RELATING Tmf PRINCIPAL: AXE 
1S*/1X,*T0 THE X, Y, and Z SYSTE'* AXES In ThAT SFQ'JF nCE */ // ) 

«ax*3 
vs 3 

*pR( 1 , nsixx 

APR ( 2 , 1 )aAP'‘ ' • . ? 

ARP(l.l)sAPR(3,nsIXZ 

ARR(?.P)alYV 

*RR(3.2)sARR(2,3)sl VZ 
A05(3,3)sIZZ 

call SY'^OL f "*AX ,n, ARR , E ,CRR , IE R^ ) 

TFdFRR ,vE, 0) GO 33? 

336 J=).3 
print 33;,j,f(..n 

537 C0P‘-aT(ix,*PIGEvvAl.F(*T1*^ s *»12,S//) 
oc!\T 339, J 

5^9 Ar5«tTfix,*EIGENVECTnt,r*T1*)*/) 

PRI .T (ARR( I , J) , Is) , 5) 

33« P'JR-iTnx,3(El«,fe,SX)////i 
536 C.r t 
GC T- 33^ 

552 PPI' T 333, TFBR 

3 53 P'-*'AT(ix,*ERRPR -• TERR s *TS) 

3 3?. jT-p 
c • "i 
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